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Abstract 

Cosmological solutions of the Brans-Dicke theory with an added cosmological 
constant are investigated with an emphasis to select a conformal frame in or- 
der to implement the scenario of a decaying cosmological constant, featuring 
an ever growing scalar field. We focus particularly on Jordan frame, the origi- 
nal frame with the nonminimal coupling, and conformally transformed Einstein 
frame without it. For the asymptotoic attractor solutions as well as the "hesita- 
tion behavior," we find that none of these conformal frames can be accepted as 
the basis of analyzing primordial nucleosynthesis. As a remedy, we propose to 
modify the prototype BD theory, by introducing a scale-invariant scalar-matter 
coupling, thus making Einstein frame acceptable. The invariacne is broken as a 
quantum anomaly effect due to the non-gravitational interactions, entailing nat- 
urally the fifth force, characterized by a finite force-range and WEP violation. 
A tentative estimate shows that the theoretical prediction is roughly consistent 
with the observational upper bounds. Further efforts to improve the experimental 
accuracy is strongly encouraged. 
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1 Introduction 



The cosmological constant is a two-step problem. First, the observational upper 
bound to A is more than 100 orders smaller than what is expected naturally from 
most of the models of unified theories. Secondly, some of the recent cosmological 
findings seem to suggest strongly that there is a lower bound as well [TJ , though 
it might be premature to draw a final conclusion. To understand the first step 
of the problem, theoretical models of a "decaying cosmological constant" have 
been proposed They are based on some versions of scalar-tensor theories of 
gravity essentially of the Brans-Dicke type M. Attempts toward the second step 
have also been made by extending the same type of theories ||. 

As a generic aspect of the scalar-tensor theories, however, one faces an inherent 
question on how one can select a physical conformal frame out of two obvious 
alternatives, conveniently called J frame (for Jordan) and E frame (for Einstein), 
respectively. The former is a conformal frame in which there is a "nonminimal 
coupling" that characterizes the Jordan-Brans-Dicke theory, but can be removed 
by a conformal transformation, sometimes called a Weyl rescaling, thus moving 
to E frame in which the gravitational part is of the standard Einstein-Hilbert 
form. 

None of realistic theories of gravity is conformally invariant. Consequently 
physics looks different from frame to frame, though physical effects in different 
conformal frames can be related to each other unambiguously. The latter fact is 
often expressed as "equivalence" 0, though sometimes resulting in confusions. 

A conformal transformation is a local change of units ||. In the context of 
Robertson- Walker cosmology, it is a time-dependent change of the choice of the 
cosmic time, measured by different clocks. In the prototype BD model with the 
scalar field decoupled from matter in the Lagrangian, the time unit in J frame is 
provided by masses of matter particles, whereas the time in E frame is measured in 
units of the gravitational constant, or the Planck mass. As will be demonstrated 
explicitly, the way the universe evolves is quite different in the two frames. We 
try to see how one can use this difference to select a particular frame. 

We confine ourselves mainly to the analysis of the primordial nucleosynthesis 
which is known to provide one of the strongest supports of the standard cosmol- 
ogy. We also focus on the simplest type of the theories which may apply only to 
the first step of the problem as stated above. The result obtained here will still 
serve as a basis of more complicated models to be applied to the second step. 

Suppose first that at the onset of the whole process of nucleosynthesis the 



2 



universe had reached already the asymptotic phase during which it evolved ac- 
cording to the "attractor" solution for the BD model with A added. We find 
that, unlike in many analyses based on the BD model without the cosmological 
constant 0, the physical result is acceptable in neither of the two frames for any 
value of u, the well-known fundamental constant of the theory. Conflicts with 
the standard picture are encountered also in the early epoch just after inflation 
and in the dust-dominated era. 

We then point out that the cosmological solution may likely show the behavior 
of "hesitation," in which the scalar field stays unchanged for some duration. If 
nucleosynthesis occurred during this phase, both frames are equally acceptable, 
giving no distinction between them. Outside the era of nucleosynthesis, however, 
we inherit the same conflicts for both conformal frames; hesitation itself offers no 
ultimate solution. 

Most of the conflicts can be avoided in E frame, as we find fortunately, if, 
contrary to the original model, the scalar field is coupled to matter in a way 
which is not only simple and attractive from a theoretical point of view but is 
also roughly consistent with observations currently available. 

In reaching this conclusion in favor of E frame, we emphasize that how a 
time unit in certain conformal frame changes with time depends crucially on how 
the scalar field enters the theory. Searching for a correct conformal frame is 
combined intricately with the search for a theoretical model which would lead to 
a reasonable overall consistency with the cosmological observations. 

In section 2, we start with defining the model first in J frame, then apply 
a conformal frame moving to E frame. We discuss in Section 3 the attractor 
solution in some detail, including elaborated comparison between the two con- 
formal frames. Section 4 discusses comparison with phenomenological aspects 
of standard cosmology, particularly primordial nucleosynthesis, dust-dominated 
universe and pre-asymptotic era. We enter discussion of the "hesitation behavior" 
in Section 5, still finding disagreement with standard cosmology. To overcome 
the difficulty we face, we propose in Section 6 a revision of the theoretical model 
by abandoning one of the premises in the original BD model, but appealing to a 
rather natural feature of scale invariance. The analysis is made both classically 
and quantum theoretically. As we find, the effect of quantum anomaly entails 
naturally the "fifth force," featuring the finite-force range and violation of the 
weak equivalence principle (WEP). Importance of further experimental studies is 
emphasized. Final Section 7 is devoted to the concluding remarks. Three Appen- 
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dices are added for some details on (A) another attractor solution, (B) mechanism 
of hesitation, and (C) loop integrals resulting in the anomaly. 



2 The model 

We start with the Lagrangian in J frame as given by 

C = sf^ (l^ 2 R ~ \tg^d^d v <\> + n A + L mattcr ) , (1) 
where our scalar field is related to BD's original notation ip = (p BD by 

<P = (2) 

also with the constant £ related to their uj by ^uj = 1/4. Notice that we use a 
unit system of c = % = 8nG(= Mp 2 ) = lfj 

We prefer <fr defined above because by doing so we write equations in a form 
more familiar with conventional relativistic field theory, avoiding the apparent 
singularity p^ 1 in the kinetic term. Also e in ([]]) can be either of ±1 or 0. 
Examples of e = — 1 are provided by the dilaton field coming from 10-dimensional 
string theory |J , and the scalar field representing the size of internal space arising 
from compactifying iV-dimensional spacetime with N > 6, while the latter with 
N = 5 gives e = 0.0 

Various models have been proposed for different choices of the function of <fi 
in the nonminimal coupling J7|. We adhere for the moment, however, to the sim- 
plest choice in (|Ij) expecting it to be applied to situations of immediate physical 
interests. The factor <p n is inserted because a cosmological constant in higher di- 
mensions may appear in 4 dimensions multiplied with some power of depending 
on the model. 

We assume RW metric, specializing to the radiation-dominated universe after 
the inflationary epoc, because it not only simplifies the calculation considerably 
but also applies to the era of nucleosynthesis. Also choosing k = we obatin the 
equations 

6^H 2 = 6^<p 2 + (p n A + p r - 6Hp, (3) 
Cp + mp = (4-n)C 2 n A, (4) 

pr+4:H Pr = 0, (5) 

-"^In units of the Planck time (2.71 x 10 _43 sec), the present age of the universe ~ 1.4 x 10 10 y 
is given by ~ 1.6 x 10 60 , while "3 minutes" is ~ 10 45 . 

2 A negative e is found to be essential in the calculations in Refs. (U carried out in J frame. 
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where 

H= h -. (6) 

a 

Now consider a conformal transformation defined by 

9*tMu = ^ 2 (x)g IMU , or dsl = tt 2 (x)ds 2 . (7) 

By choosing 

n = ^/ 2 0, (8) 

we rewrite (Tl|) as the Lagrangian in E frame: 

£ = V^F* (^R* - ^f<V^a + v{a) + L* matter ) , (9) 

where 

= T 1/2 e CCT , (10) 
with the coefficient ( as defined by 

r 2 = 6 + er 1 = 2(3 + 2e^), (11) 

and 

V(a) = Q- 4 f l A = A^- n / 2 e (n - 4)?,J . (12) 

Notice that the A term in ([!]) now acts as a potential even with n = 0. 

We point out that the canonical field a is a normal field (not a ghost) if the 
right-hand side of ( |TT]) is positive: 

6 + e^~ 1 >0, (13) 

even with e = — 1. This positivity condition is satisfied trivially for any finite 
value of £ if e = . The condition (|l^) should be obeyed even in the analysis in J 
frame; in the presence of the nonminimal coupling that causes mixing between 
and the spinless part of the tensor field, the sign of the total scalar-tensor sector is 
not determined solely by e. We notice that the conformal transformation serves 
also as the relevant diagonalization procedure. We impose £ > 0, to keep the 
energy of the tensor gravitational field positive. 

With the RW metric also in E frame, the cosmo logical equations are 

3if 2 = p* CT + p„, (14) 
a + 3H*a + V'(a) = 0, (15) 
p\ r + AH*p* r = 0, (16) 
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where 

P*„ = \v 2 + V{a). (17) 

An overdot implies a differentiation with respect to t*, the cosmic time in E frame 
defined by [§ 

dU = tldt, (18) 
a* = fla, (19) 

where the scale factor a* has also been introduced to define H* = a* /a*. 



3 Attractor solution 

The solutions of the cosmological equations can be obtained more easily in E 
frame than in J frame. We in fact find a special solution; 

a*(t*) = t l J\ (20) 
2 

a(Q = a + - C mt*, (21) 

4 — n 



P* r (t*) = P*rot* , with p 
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*r0 
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where the constant a is given by 



4-n 



2 

c 2 



(22) 



Ae( "" 4)CS = 7 7^^ /2 C" 2 - (23) 

(ra — A) 1 

We also find 

M**) = 7-A^C-\- 2 , (24) 
(n — A) z 

which is independent of A. 

It is easy to see that, since V(a) is an exponentially decreasing function if 
n < 4, a is pushed toward infinity, hence giving p w ~ t~ 2 showing that A e g = p w 
does decay with time even with n = 0, a purely constant cosmological constant 
in 4 dimensions |§ . According to (|2l|) the scalar field continues to grow; a unique 
feature that distinguishes our model from other models |7j] with the scalar field 
designed to settle eventually to a constant .0 

To ensure a natural condition p„. > we impose 

< 2 >J4 ± W < 25) 

3 If A is added to this type of models, this constant would yield a too large A e ff = constant 
which we attempt to avoid without fine-tuning parameters. 
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Figure 1: An example of the solutions to (|I4|)-(|T6D with n = 0. We chose A = 
( = 1.0. The initial values at logt*i = 10 are a = 12.0, da /dint* = 1.0, t^p^ = 
0.2. All the curves tend to the asymptotic lines given by fl2"0|)-(|2"2"D quickly after 
logt* ~ 20. The present age of the universe is ~ 10 60 . 



Though there is a subtlety as will be discussed in Appendix A, this is one of the 
peculiar outcomes of including A. According to (|ll|) this predicts uj too small to 
be consistent with the currently accepted lower bound. Later, however, we will 
come to propose a revised theoretical model in which ( is not related directly to 
phenomenological parameters. We also notice from (|TT|) that £ > 1/ v^6 can be 
realized only for e = — 1.0 

It is already known in the literature that the solution (|20|) - (|24"D is an attractor, 
though many of the investigations have been aimed at the question if it leads to 
a sufficient inflation |9-15|. In contrast we are interested in how the cosmological 
constant is relaxed in a manner consistent with other known aspects of cosmol- 
ogy^ For this reason we start integrating the equations sometime after the end 
of inflation when sufficient amount of matter energy is created due to reheating, 
though its details have not been well understood. Fig. 1 shows how a typical 



solution of (|i4|)- (|l6|) (with n = 0) tends asymptotically to the attractor behavior 
for logi*,>20. We emphasize that the behaviors are essentially the same for any 
values of n as far as n < 4; according to ( |i~2"D different n is absorbed into different 

c 

In this example we chose the initial value of a at the time £*i(= 10 10 ) in 

4 This is equivalent to choosing a negative uj according to some authors |fj^Jl5f . 
5 See also Ref. Ml for some of our earlier results. 
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such a way that the resultant solution varies "smoothly" around t*\. To be more 
specific, we require t^H*, the local effective exponent in a*(i*), to remain of the 
order unity. This is in accordance with assuming that the "true" initial condition 
on the more fundamental level is given at a much earlier time. 

The same attractor solutions can be obtained also in J frame. This can be done 
most easily by substituting (p0|)- (p4|) into ( |IBD and ([19]). After straightforward 
calculations we find 



A ^(2-n)/(4-n) ; tf 

A 2 lnt*, if n = 2, 



(26) 



and hence 



t a , with a = if n ^ 2 , 

ait) = { (27) 

exp (-A 2 r 1/2 a) , if n = 2, 



where A n is a constant given in terms of A, £ and n. 

Of special interest is the choice n = 0, a purely constant cosmological constant 
in 4 dimensions, giving 

i = A t l J 2 , (28) 
cs(t) = const. (29) 



The asymptotic behavior fl28|) is seen in Fig. 2 for the same solution as in Fig. 
1, under the condition t\ = t*i at the "initial" time. Notice that this behavior is 
preceded by a period of an extremely slow take-off of t(t*), which can be traced 
back to a large value of Q appearing in (|18D , related to the remark stated above 
on the initial value er(£i*). Fig. 3 shows how the asymptotic behavior ( ^9"D is 
reached. A closer look at the curve toward the initial time reveals an oscillatory 
behavior too rapid to be drawn here exactly. This might have a disturbing effect, 
as will be discussed later. 

Also for n = 2, for which some of the exact solutions have been obtained 
|fll,|14,T5i, we find an exponentially contracting universe, though this can be 
converted to an exponentially expanding universe by reversing the direction of 
time. 



The exponent a in the first of (|27|) is larger than 1 for n > 2, while it is 
negative for < n < 2. The exponent never reaches 1/2 for any finite value of n; 
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Figure 2: The relation between two times, t* in E frame and t in J frame, com- 
puted for the example shown in Fig. 1. 

getting 0.45, for example, requires n — — 18. We also point out that the behavior 
differs considerably depending on whether the "ordinary" matter (p r or p„) is 
included or not. 

Summarizing, we expect substantial differences in the behavior of the scale 
factor not only between the conformal frames but also among different values of n 
in J frame. This demonstrates how crucial it is to select a right conformal frame 
to discuss any of the physical effects. In the next section we analyze the physical 
implications to be compared with the results of standard cosmology. 

4 Comparison with standard cosmology 
4.1 Nucleosynthesis 

According to the standard scenario, light elements were created through nuclear 
reactions in the radiation-dominated universe with the temperature dropping 
proportional to the inverse square root of the cosmic time. The whole process is 
analyzed in terms of nonrelativistic quantum mechanics in which particle masses 
are taken obviously as constant. In scalar-tensor theories, on the other hand, par- 
ticle masses depend generically on the scalar field, hence on time. The prototype 
BD theory is unique in that masses are true constants thanks to the assumption 
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Figure 3: tH in J frame, computed for the example shown in Fig. 1, showing 
that the scale factor tends to a constant. Exact behaviors are omitted for t very 
close to the initial time. See Fig. 4 for more details. 

that the scalar field is decoupled from matter in the Lagrangian. Q For this reason 
the mass m of a particle is a pure constant in J frame, specifically denoted by 
m Q . 

We should recall that we have so far considered relativistic matter alone. One 
may nevertheless include massive particles which play no role to determine the 
overall cosmological evolution, but may serve to provide standards of time and 
length. 

To find how masses become time-dependent after a conformal transformation, 
it is sufficient to consider a toy model of a free massive real scalar field $, described 
by the matter Lagrangian in J frame; 

Matter = (-^<9^$ - ^71^ , (30) 

where no coupling to is introduced. After the conformal transformation (0), we 
find 

^matter = (-^Q^d^d^ - . (31) 

The kinetic term can be made canonical in terms of a new field defined by 

$ = (32) 

6 This is a major point that distinguishes between BD model and the one due to Jordan p"7| ] 
who pioneered the nonminimal coupling. 
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thus putting ( pip into 

1 



where 



^matter = ( --g<?V^V v $* - -m& ) , (33) 



Z>„$. = [<^ + C(<V)]$*, (34) 

m 2 = fi~ 2 m 2. (35) 

We point out that this relation holds true generally beyond the simplified model 
considered above. 

From (D, (|10]) and © we find 

hence giving a time- dependent mass in E frame. This would result in the reduction 
of masses as much as 1 — 1/a/10 ~70% if n = 0, for example, in the period 100- 
1000 sec, during which major part of synthesis of light elements is supposed to 
have taken place with the temperature dropping in the same rate. Obviously this 
is totally in conflict with the success of the standard scenario. 

In this way we come to a dilemma; J frame is selected uniquely because of 
constancy of masses as taken for granted in conventional quantum mechanics to 
analyze the physical processes, while E frame is definitely preferred to have the 
universe that cooled down sufficiently for light elements to form. 

In passing we offer a simple intuitive interpretation on the relation between 
the two conformal frames. In the present context, the time unit r in E frame 
is provided by m^ 1 ~ Q. The time t measured in units of r may be defined by 
dt = dt*/T ~ Q~ 1 dt*. Comparing this with (18[), we find t = t. On the other 



hand, the only dimensionful constant in E frame is Mp. In this sense E frame 
corresponds to the time unit provided by the gravitational constant. For n = 0, 
for example, the microscopic length scale provided by m~ l expands in the same 
rate as the scale factor a*(i*), hence showing no expansion in a(t) in dl9|) . 

4.2 Dust-dominated era 

In the standard theory the radiation-dominated era is followed by the dust- 
dominated universe. Its description is, however, likely problematic, as will be 
shown. 
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Due to the equation 

□*<7 = cCT*, (37) 
with T* = — p^, the right-hand side of (|T^) acquires an additional term 

c(p*d, (38) 

where, for the non-relativistic matter density p^, c = 1 for the prototype BD 
model, but we allow c to be different in the proposed revision which will be 
discussed later. Suppose the total matter density is the sum p* r + p* d , which 
would replace p* r in ([14]). Corresponding to (0), we find 

p\ d + 3#*p* d = -cCp* d (T, (39) 

where the right-hand side is included to meet the condition from the Bianchi 
identity 

We obtain the attractor solution with 

4 - c 

aJU)=t°', with a*= — — , (40) 

o 

and 

4 — c 1 9 j . 

P*do = — - ~ A C 2 - (41) 

Equation fl2~l~|) still holds true. Notice, however, that (|40"D gives a* = 1/2 for 
c = 1, the same behavior as in the radiation-dominated era. This seems "uncom- 
fortable," if we wish to stay close to the realm of the standard scenario, though 
we may not entirely rule out a highly contrived way for a reconciliation. On the 

other hand, we would obtain the conventional result a* = 2/3 for c = 0. 

i li 

We also find that the relation t ~ t* remains unchanged (for n — 0), and 
then follows a = const in J frame again for c = 1. This is certainly disfavored, as 
in the analysis of nucleosynthesis. 

4.3 Pre-asymptotic era 

As we noticed in Fig. 3, H in J frame seems to be oscillatory around zero in early 
epochs, arousing suspicion of a contracting universe. More details towrad t\ = t*\ 
are shown in Fig. 4, in which tH is plotted against logt* instead of logt; the dip 
of tH would be too sharp to be shown if plotted against logt. We also plotted 
log a, which does show a decrease of the scale factor in J frame. The example here 
shows that the scale factor a which comes to at rest toward the asymptotic era is 
even smaller than at £*i. The exact amount of contraction depends on the choice 
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of the parameters, still making it considerably difficult to reach a compromise 
with the idea that the early universe had cooled down sufficiently to trigger the 
process of nucleosynthesis. 
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Figure 4: More detailed behavior of tH and log a plotted against logi*, computed 
for the example shown in Fig. 1. The value of the scale factor a in J frame when 
it comes to at rest toward the asymptotic behavior is even smaller than a± at the 
"initial" time. 

5 Hesitation behavior 

We have so far concentrated on the attractor solution, to which some solutions, 
depending on the initial conditions, do tend smoothly, as demonstrated in Fig. 1, 
for example. We point out, however, there is an important pattern of deviation 
from this smooth behavior. As illustrated in Fig. 5, the scalar field may remain 
almost at rest temporarily before entering the asymptotic phase in which it re- 
sumes to increase to approach the attractor solution. This "hesitation" behavior 
may occur if p* r <C p w at the initial time, as will be shown in detail in Appendix 



Notice that t^H*, which equals the effective exponent a* if the scale factor is 
approximated locally by a*(t*) ~ tends to 1/2 after some wiggle-like behavior 
that separates the plateau of the same value 1/2 in the hesitation period. See 
also Appendix B for the mechanism behind another plateau of = 1/3. 

With the scalar field nearly constant during this hesitation phase, particle 
masses are also nearly constant, and all the other cosmological effects are virtually 
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Figure 5: An example of a hesitation behavior, obtained by choosing the initial 
value of £*p* r as 2.0 x 10~ 15 , with other parameters the same as in the solution 
in Fig. 1. The i^-domination occurs for 11 <;logt* <; 30, while the hesitation 
behavior is seen clearly for 35 ^logi* ^49, followed by the usual radiation- 
dominated universe. Notice that "3 minutes" corresponds to ~ 10 45 . 

the same as those in the standard theory as long as p* CT <C p* r , as is the case for 
logt*,>35 in the example of Fig. 5, in which we chose the parameters in such a 
way that the hesitation period logt* = 35-48 covers the era of nucleosynthesis. 
Moreover, the constant scalar field makes the conformal transformation fllBp and 
( |T5| ) trivial, implying that the two conformal frames are essentially equivalent to 
each other. 

In more details, however, we find some differences between them. In our 
example, we carried out the transformation ( |T8"D and (0), showing first in Fig. 
6, how the two time variables £* and t are related to each other; the period of 
t ~ i* during hesitation is present in addition to the behavior as shown in Fig. 2 
without the hesitation behavior. 

As in Fig. 4, we find in Fig. 7 that the universe in J frame had experienced 
a considerable contraction prior to the epoch of nucleosynthesis, making the sce- 
nario of the evolution in early epochs desperately different from the standard 
theory. See Appendix B for its origin. 

No such problems will occur if we are still during the hesitation period at the 
present time. If this happens, however, no distinction is present between the two 
frames, offering no obstruction to choose E frame either. From this point of view, 
we do not consider this possibility any further. 
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Figure 6: The relation between the two times for the solution shown in Fig. 5. 

1/2 

The behaviors t ~ and t ~ t* are seen for 35 <^logt* ^49 and logt* ;>50, 
respectively. 



One might suspect that all of these "conflicts" with the standard picutre come 
directly from the theoretical models to start with. In fact it seems obvious that 
we would be in a much better position if particle masses were constant in E frame 
rather than in J frame. This can be achieved, as we will show, by modifying one 
of the assumptions in the prototype model in a rather natural manner. 



6 A proposed revision 
6.1 Classical thoery 



Unlike in the original BD model, let us start with the matter Lagrangian in J 
frame: 

Matter = (-~<T W$ - \f<?& - ^A$ 4 ) , (42) 

in place of fl30]). We introduced the coupling of to the matter field $ 
hence abandoning one of the premises in the original BD theory. Also there is no 
mass term of $; "mass" of the field <3> is f(f) which is no longer constant. With the 
choice (^), therefore, J frame loses its privilege to be a basis of the theoretical 
analysis of nucleosynthesis. We also introduced the self-coupling of $ to illustrate 
the effect of the quantum anomaly. 

After the conformal transformation we obtain the same Lagrangian (|33"D (plus 
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Figure 7: Like in Fig. 4 for the solution without hesitation, the solution shown 
in Fig. 5 results in the considerable amount of contraction of the scale factor in 
J frame in the early universe. 



the self-coupling term) but with replaced by 

m% = fVV- 2 = r 1 / 2 , (43) 

which is obviously constant, making E frame now a relevant frame for realistic 
cosmology. We point out that no matter coupling of a is present in E frame. 
Absence of the coupling in the Lagrangian having no nonminimal coupling im- 
plies a complete decoupling, unlike the corresponding situation in J frame in 
the prototype model. This corresponds to the choice c = in (|37|) - (f|l|), thus 



leaving the results in E frame the same as those of standard cosmology also in 
dust-dominated universe. On the other hand, the scalar field can be detected in 
no ways by measuring its contribution to the force between matter objects, or 
the conventional tests of General Relativity, hence removing the constraints on 
uj (or £) obtained so far. Its effect may still be manifest through cosmological 
phenomena. 

The time scale in E frame is provided commonly by particle masses and the 
gravitational coupling constant. This implies that no time variability of the grav- 
itational constant should be observed if measured by atomic clocks with their 
unit given basically by particle massesj] 

7 Strictly speaking, the time unit of atomic clocks depends also on the fine-structure constant, 
whose constantcy is assumed, however, at this moment. 
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The above scheme is attractive because the coupling constant / is dimension- 
less, hence vesting scale invariance in the gravity-matter system except for the A 
term. By applying Noether's procedure we obtain the dilatation current as given 
by 

2(-%ct + (d u + 2(d„<j) $2 



(44) 



which is shown to be conserved by using the field equations. This conservation 
law remains true even after the conformal transformation, with a nonzero mass 
as given by ((43[). This implies that the scale invariance is broken spontaneously 
due to the trick by which a dimensionful constant Mp(= 1) has been "smuggled" 
in (||). In this context a is a Nambu-Goldstone boson, a dilaton. This invariance 
together with constancy of particle masses will be lost, however, if one includes 
quantum effects due to the non-gravitational couplings among matter fields, as 
will be sketched below. 





(a) 



(b) 



Figure 8: One-loop diagrams giving an anomaly. The blobs are for the coupling 
(A/4!)$ 4 , while the crosses represent a coupling to a, the linear term in 



6.2 Quantum anomaly 

Consider one-loop diagrams for the coupling between $^ and a assumed to carry 
no momentum, as illustrated in Fig. 8, arising from the non-gravitationa cou- 
pling (A/4!)$ 4 . They will result, according to quantum field theory, generally in 
divergent integrals, which may be regularized by means of continuous spacetime 
dimensions. Corresponding to this, we rewrite previous results extended to N 
dimensions. Equation (|32l ) is then modified to 

$ = jy- 1 ^, (45) 

where v = N/2. Also the relation ([8]) is changed to 

n^V^-VA"-!). (46) 
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The Lagrangian ([42]) is now put into 

^matter = [—g^V^V^ ~ S^j , (47) 

where T>^ is given by (|34j ) while 5* is defined by 

3 = e 2{ »- 2Ka S , (48) 

with 

5 = r X y^ + ^Ae (49) 

As a result the right-hand side of (|43|) is multiplied by e^ -4 -*^ 7 , which might 
be expanded with respect to (N — 4)(a according to ([[(]); decoupling occurs only 
in 4 dimensions. Following the rule of dimensional regularization, we keep N off 
4 until the end of the calculation including the evaluation of loop integrals which 
are finite for N ^ 4. 

The divergences coming from the loops are represented by poles (N — 4) _1 
which cancel the factor N — A that multiplies the a coupling as stated above, 
hence yielding an effective interaction of the form 

" LU = 9^>, (50) 
for iV — > 4, where the coupling constant g a is given by 

g* = C^a, (51) 

with the details of computation presented in Appendix C. In ([j0|) we reinstalled 
Mp 1 = y/8nG to remind that this coupling is basically as weak as the gravita- 
tional interaction. 

Suppose field at rest is a representative of dust matter. We may then take 
p* d ps ml §l. Comparing fl50|) with fl37|) we find 

g <y = ~c(. (52) 

In this way we obtain g a , and hence c, which are nonzero finite due to a non- 
gravitational interaction. 

It should be warned that our determination of c in (|37|) never implies that a 
couples to the trace of the energy-momentum tensor. From (^) we find that the 
coefficient g a depends on A which is not related to the mass directly, hence may 
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differ from particle to particle. If we include many matter particles, the right- 
hand side of (jT?|) should be the sum of corresponding components with different 
coefficients. For this reason the force mediated by a fails generically to respect 
WEP. This was shown more explicitly in our QED version [|T9[|. 

We notice that the above calculation is essentially the same as those by which 
various "anomalies" are derived, particularly the trace anomaly PC| . The rele- 
vance to the latter can be shown explicitly if we consider the source of a in the 
limit of weak gravity. 

We first derive 

<V M = 2C(^ - 2)V=^s; (53) 

indicating that scale invariance is broken explicitly for iV 7^ 4. Now deriving (|50|) 
is essentially calculating the quantum theoretical expectation value of So between 
two 1-particle states of $*; 

<%(v- 2)S >*.= g a ml , (54) 

ignoring terms higher order in £(= (Mp 1 ). Combining this with (|37D we may 
write 

a** =< <v M >** • ( 55 ) 

This equation shows that a couples to breaking of scale invariance effected by 
quantum anomaly,^ though this simple realtion is justified only up to the lowest 
order in Mp . 

One might be tempted to extend the analysis by identifying $ with the Higgs 
boson in the standard electro-weak theory or Grand Unified Theories, hence pre- 
dicting observable consequences. As we find, however, the realistic analysis is 
ought to be more complicated; we must take contributions from other couplings 
including the Yukawa and QCD interactions into account. Even more serious is 
that we are still short of a complete understanding of the "content" of nucleons, 
the dominant constituent of the real world. We neverthelss attempt an analysis, 
as will be sketched briefly, leaving further details to the future publications. 

From a practical point of view, we need the coupling strength of a to nucleons, 
through the couplings to quarks and gluons. By a calculation parallel to that 

8 In this respect we rediscover the proposal by Peccei, Sola and Wetterich pH , though their 
approach to the cosmological constant problem is different from ours. 
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leading to fl5"T|) we obtain^ 

c = -™L ~ 0.3, (56) 

7T 

for a ex-quark coupling, where a s is the QCD analog of the fine-structure constant, 
most likely of the order of ~ 0.2. We then evaluate 

c N m N = c q <^2miq*iq*i >n, (57) 



where the nucleon matrix element has been estimated to be ~ 60MeV, |22 which 
together with (|56| ) would give 

c N ~ 1.8 x 10~ 2 . (58) 

It is interesting to note that this is rather close to the constraints obtained from 
observations, as will be shwon. 

6.3 Fifth force and cosmology 



By replacing £ in ( |lT|) by c£ also choosing e = +1 in accordance with the con- 
ventional analysis, we obtain 

C - 2 C 2 = 8(3 + 2c;), (59) 

or 

< ~ <0.8 x 10~ 2 , (60) 



if we accept tu^lO 3 obtained from the solar-system experiments, assuming the 
force-range of a longer than 1 solar unit. Notice that a is now a pseudo Nambu- 
Goldstone boson which likely acquires a nonzero mass. Combining this with the 



constraint (25), we find 

\c\ <1.6 x 10" 2 , (61) 

which we find is nearly the same as (|51f). 

Similar constraints may come from the "fifth force" phenomena which are 
characterized by a finite force-range and WEP violation both of which are generic 
in the present model. The parameter a^, the relative strength of the fifth force, 
is given by 

« 5 ~ 2g 2 a ~ 4, (62) 

9 The coefficient 6 in eq. (52) of Ref. should be replaced by 15/4. (Also (3 is our present 
C-) The QCD result, our (p6[), is obtained by multiplying further by the factor (TV 2 — 1)/2N 
which is 4/3 for N = 3. 
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expecting that the coupling comes mainly from the one to nucleons. 

Since g a may depend on the object a couples to, as was pointed out, a§ may 
also depend on the species of nuclei, for example, between which a is exchanged, 
to be denoted by a 5i j. 

From the observations carried out so far, we have the upper bounds given 
roughly by [f23 

|a 5ij |<i(r 5 . (63) 

In view of the fact that the result depends crucially on the assumed value of the 
force-range as well as the model of WEP violation, we may consider that the 
estimate (|62|) with (|58|) is approximately consistent with (|63"D , hence providing a 
renewed motivation for further studies of the subject both from theoretical and 
experimental sides. 

The analysis is still tentative particularly because (|56"D is justified only to the 
lowest order with respect to a s , though the renormalization-group technique can 
be used to include the leading-order terms. Potentially more important would be 
to estimate the contribution from gluons. Corresponding to the second term of 
(|C.1|) , we should include the direct coupling of a to the QCD coupling constant 
g s as given by 

lim ((u - 2)g s Z g a = 2ir(f3 g s a s a, (64) 

where (3q = (47r)~ 2 (ll — 2rif/3) with nf the number of flavors. It is yet to be 
studied how this coupling would affect the simple result (|58|) through the gluon 
content of a nucleon. 

It should also be emphasized that the right-hand side of fl57|) (even with the 
modification stated above) is quite different from the matrix element of the con- 
ventional energy-momentum tensor or its trace related directly to observations; 
only the anomalous part participates.^ 

One might argue that cjy which would be too large to be allowed by the phe- 
nomenological constraints could emerge if we apply the same type of calculation 
to a nucleon considered to be an elementary particle, as was attempted in the 



simpler QED version [19]. We point out, however, the finite size of a composite 
nucleon would serve to suppress ultra-violet divergences, thus failing to produce 
an anomaly, which is a manifestation that the underlying theory is divergent. 

It is rather likely that a couples also to the nuclear binding energy which is 
generated supposedly by the exchanged mesons. This would make the analysis 



10 An analysis due to Ji [Q showing that the trace anomaly contributes approximately 20% 
of the nucleon mass might also be suggestive. 
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of composition-dependence even more complicated .0 

We now turn to cosmological aspects. Adding fl5"0|) to the "classical" mass 
term, we typically obtain 

- L' m t = -m*$l, with ml = ml (l-c(a + ■■■). (65) 

We focus on the cosmological background <x(£*) rather than the sapce-time fluc- 
tuating part which would mediate a force between matter objects as considered 
above. In this sense m* depends on time. We must then apply another confor- 
mal transformation to cancel this effect. If, however, c is sufficiently small, as 
indicated in (|6TD, one may expect that the expansion in (|65"D is exponentiated 
giving 

m*(cr) « m, e- cC,T ~ C /2 , (66) 

where we have used the asymptotic behavior fl2"T| ) with n — 0. We expect fl5"E|) 
holds true approximately for realistic nucleons or nuclei. 

Notice that the final expression on the time variation is independent of £. 
The resulting conformal frame is expected to be close to E frame. A small |c| is 
also favored from (HO) for the attractor solution in the dust-dominated unvierse. 



It would further follow that G/G is somewhat below the level of 10 -10 y _1 , in 
accordance with the observations [^(| . 

If, on the contrary, |c| is "large," we may even not be able to compute the 
required conformal transformation unless we determine higher-order terms in the 
parenthesis in (p5]). All in all, we would be certainly "comfortable" if |c| is 
sufficiently small. On the other hand, it seems unlikely that |c| is smaller than 
unity by many orders of magnitude because we know no basic reason why it should 
be so. The present constraint fl6"Tj ) might be already close to the limit which one 
can tolerate in any reasonable theoretical calculation. In this sense, probing 
a 5 with accuracy improved by a few orders of magnitude would be crucailly 
important to test the proposed model of broken scale invariance. If we come to 
discover any effect of this kind, it would provide us with valuable clues on how 
nucleons and nuclei are composed of quarks and gluons. 



11 Most of the past analyses on the composition-dependent experiments |25[| have been made 
based on the assumption that the fifth force is decoupled from the nuclear binding energy, which 
might be too simplified from our point of view. 
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7 Concluding remarks 



Having introduced a scalar field in order to relax the cosmological constant within 
the realm of the standard scenario, we come to a conclusion: At the classical level, 
E frame is the only choice provided the J frame version has a scale invariant cou- 
pling between the scalar field and the matter fields without the intrinsic mass 
terms. Most crucial are constancy of particle masses and the expansion law of 
the universe during the epoch of primordial nucleosynthesis. A quantum anomaly 
serves naturally to break scale invariance explicitly, offering yet another support 
for the occurrence of the fifth force featuring WEP violation. A tentative cal- 
culation based on QCD yields the coupling strength roughly consistent with the 
observational upper bounds. The physical conformal frame should then remain 
close to E frame. Improved efforts to probe the fifth force is encouraged, though 
detailed theoretical predictions are yet to be attempted. 

We point out, however, there is a possible way to leave a completely decoupled 
even with quantum effects included. We may demand that the 0-matter coupling 
in (^2) has a coupling constant / which is dimensionless in any dimensions. This 



can be met if we replace the second term in the parenthesis in ( [42]) by 

- \f^ v - X) ^\ (67) 

resulting in m* which is shown to be completely cr-independent in E frame, even 
with the quantum effect included. 

We find, on the other hand, that the term of the nonminimal coupling of the 
form (p 2 R, as in the first term of ([!]), is multiplied always with a dimensionless 
constant for any dimensions. This is a fact that underlies the whole discussion of 
scale invariance, making another difference from the models |7| which allow more 
general functions of the scalar field. Scale invariance is respected also in a new 
approach to the scalar-tensor theory based on M4 x Z2 |27J . 

The present study is limited because the model with a single scalar field might 
be too simple to account for a possible nonzero cosmological constant. Q Looking 
further into the hesitation behavior would be still useful to acquire more insight 
into the time- (non) variability of various coupling constants, probably a related 
issue which seems to deserve further scrutiny [ITH]. 



12 See, for example, || for a model with another scalar field implementing a "sporadic" oc- 
currence of a small but nonzero cosmological constant. 
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Appendices 
A Another attractor solution 

It seems interesting to see what the fate of the matter energy p„ r would be if 
( violates (|25|). Rather unexpectedly, p„ < is evaded automatically. The 
solutions for ( smaller than the critical value given by (|25|) tend to the "vacuum" 
solution given by (n = for simplicity) 

a(t,) = t a * with a, = ^C" 2 , (A.l) 

o 

a (t») = — hiU + a, (A.2) 

P*r(*0 = °> ( A -3) 



64C 4 



where a is given by 



Notice that these agree with (^Oj)-(H) with n = for ( = £ = 1/2. 

The two solutions may be depicted schematically as in Fig. 9. For ( > ( , 
the solution with p* r > 0, represented by Track 1, is an attractor. Suppose we 
descend along Track 1. At ( = £ , one switches to Track 2 for p„, = instead 
of yielding negative matter energy. We in fact have examples to show that the 
lower half of Track 1 is a repeller. 

We may respect (^5l) as far as we are interested only in the universe which 
accommodates nontrivial matter content asymptotically. 
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Figure 9: The two Tracks 1 and 2 represented by two straight lines crossing each 
other correspond to the asymptotic solutions (p2|) and ( |A.3| ), respectively. The 
horizontal axis is (~ 2 , while the vertical axis stands for p* r o. The upper half of 
Track 1 and the right half of Track 2 represent attractors, other portions repellers. 
A "switching" occurs at the "crossing" at £5~ 2 - 

B Hesitation behavior 

Suppose there was a period in which p* a ^> p* r in the very early universe. This 
is reasonbaly expected if reheating after inflation was not too much sufficient to 
recover p* r which had been extremely red-shifted during inflation. Notice that 
p* CT should have stayed basically of the order of A when a was rolling down the 
slope of the exponential potential as given by ([12]). 

Also the exponential slope is so steep that V(a) becomes small rapidly, leading 
to the -fT^-dominated unvierse, where 



As is well known only the kinetic energy of a scalar field as the matter density is 
equivalent to the equation of state p — p, resulting in the expansion law 



(B.l) 




(B.2) 



Using this in (|15 ) also with V ignored we find a solution 



o ~ Pi' 1 , 




and hence 



a ~ /31nt*. 



(B.4) 
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It also follows 

P* CT « K a ~ i/3 2 t; 2 . (B.5) 

The coefficent /3 can be determined if we substitute ( |B.2|) and (BJ5) into ( |l4|) also 
with V(cr) dropped: 

/3=/|. (B.6) 

In fact the solution given by (|B.2|) and (|B.4|) with flB.6j ) is found to be an attractor. 
On the other hand, substituting ( p.2| ) in (|16"D yields 

P*r ~ ^ 4/3 , (B.7) 

which falls off more slowly than p^ as given by ( |B.5|) . For this reason p* T soon 
catches up and eventually surpasses p* CT . After this time, the universe enters the 
radiation-dominated epoch with 

H* = (l/2)t~\ (B.8) 

Using this in (|15|) we obtain 

& ~ t; 3 / 2 , (B.9) 

showing that cr comes to at rest quickly. This is the beginning of the hesitation 
phase. We also find 

p*,« K a ~ t ; 3 , (B.10) 

which decreases much faster than p„. ~ t^ 2 . 

The potential V(cr) ~ e -4 ^ 7 had already been very small at the onset of the 
hesitation period, staying there since. Eventually, however, K a reaches this small 
value, so that ( |T5| ) has to be solved with V included again, hence the end of the 
hesitation. 

Fig. 5 shows a period in which ( |B.2|) is obeyed. We can also derive the 
behavior of a in J frame during the ^-dominated epoch. From ( |B.4| ) we first 
find 

(B.ll) 

Using this together with (|B.2j) in (|19D we obtain 



a(Q~t & , with a=--C/3. (B.12) 
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From (|25T) with n = and (]B.6f) we find 



a < - - -p = -0.108 < 0, (B.13) 
3 V6 

implying contraction of the universe in J frame during the -K^-dominated universe. 
In Fig. 7 we recognize the slope of log a against logi* for Q = 1; 

a = I - yf = -0.483. (B.14) 

C Loop integrals 

We show some of the details of calculating loop diagrams of Fig. 8. The basic 
coupling to the linear a is obtained by expanding the exponential in ([48]): 



- LU = 2((u - 2)a Qm^ + , (C.l) 

where m^ = ^ _1 / 2 . The diagram (a) comes from the second term of ( |C1| ) while 
the diagram (b) from the first. The 4-vertex in (b) represents the simple (A/4!)$^ 
without a attached. 

Considering a simply as a constant, we compute the amplitude for the diagram 

(a): 

M a = ~t(27r)- N 2(u - 2)(o-\ I d N k 1 2 (C.2) 

We use 

dN \ k 2 + m 2 o) = ^V r(i - u), (c.3) 



where we have put N = 4 except inside the T function. Substituting ( |C.3| ) into 
( |C.2| ), and further using 

(2 - u)T(l - v) = -^—{2-v)T{2-v) 

' r(3-z/)^-l, (C.4) 



1 - v 

we obtain 



M a = (a-^ml . (C.5) 
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Corresponding to (b), we find 



M b = z(2n)- N 2(u - 2)(aml Q \ / d N k — - (C.6) 
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The integral is now obtained by operating (— <9/<9m^ ) to the integral in ( |C.3| ), 
thus yielding 

M b = M a . (C.7) 
Adding M. a and A4b, we finally obtain 

M = (*^- 2 m%, (C.8) 
from which follow (|50|) and (51). 
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